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Abstract. In this article, we review the construction and properties of 
some popular approaches to modeling LIBOR rates. We discuss the fol- 
lowing frameworks: classical LIBOR market models, forward price mod- 
els and Markov-functional models. We close with the recently developed 
afhne LIBOR models. 

1. Introduction 

Interest rate markets are a large and important part of global financial 
markets. The figures published by the Bank for International Settlements 
(BIS) show that interest rate derivatives represent more than 60% of the 
over-the-counter markets over the years, in terms of notional amount; cf. 
Table 1. Hence, it is important to have models that can adequately describe 
the dynamics and mechanics of interest rates. 

There is a notable difference between interest rate markets and stock or 
foreign exchange (FX) markets. While in the latter there is a single under- 
lying to be modeled, the stock price or the FX rate, in interest rate markets 
there is a whole family of underlyings to be modeled, indexed by the time 
of maturity. This poses unique challenges for researchers in mathematical 
finance and has led to some fascinating developments. 

The initial approaches to interest rate modeling considered short rates 
or instantaneous forward rates as modeling objects, and then deduced from 
them tradable rates. More recently, effective rates, i.e. tradable market rates 
such as the LIBOR or swap rate, were modeled directly. Models for effec- 
tive rates consider only a discrete set of maturity dates, the so-called tenor 
structure, which consists of the dates when these rates are fixed. A review 
of the different approaches to modeling interest rates is beyond the scope 
of the present article. There are many excellent books available, focusing 
on the theoretical and practical aspects of interest rate theory. We refer the 
reader e.g. to Bjork (2004), Musiela and Rutkowski (1997), Filipovic (2009), 
or Brigo and Mercurio (2006). 

The aim of this article is to review the construction and basic properties 
of models for LIBOR rates. We consider the following popular approaches: 
LIBOR market models, forward price models and Markov-functional models, 
as well as the recently developed class of affine LIBOR models. In section 3 
we will present and discuss some basic requirements that models for LIBOR 
rates should satisfy. These are briefly: positivity of LIBOR rates, arbitrage 
freeness and analytical tractability. 
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Dec 2006 Dec 2007 Jun 2008 Dec 2008 



Foreign exchange 



40,271 56,238 62,983 49,753 

291,581 393,138 458,304 418,678 

7,488 8,469 10,177 6,494 

7,115 8,455 13,229 4,427 

28,650 57,894 57,325 41,868 

43,026 71,146 81,708 70,742 



Interest rate 

Equity-linked 

Commodity 



Credit default swaps 
Unallocated 



Total 



418,131 595,341 683,726 591,963 



Table 1. Amounts outstanding of over-the-counter (OTC) 
derivatives by risk category and instrument (in billions of US 
dollars). Source: BIS Quarterly Review, September 2009. 



There are two natural starting points for modeling LIBOR rates: the rate 
itself and the forward price. Although they differ only by an additive and a 
multiplicative constant, cf. (2.1), the model dynamics are noticeably differ- 
ent, depending on whether the model is based on the LIBOR or the forward 
price. In addition, the consequences from the point of view of econometrics 
are also significant. 

Modeling LIBOR rates directly, leads to positive rates and arbitrage- 
free dynamics, but the model is not analytically tractable. On the other 
hand, models for the forward price are analytically tractable, but LIBOR 
rates can become negative. The only models that can respect all properties 
simultaneously are Markov-functional models and afRne LIBOR models. 

The article is organized as follows: in section 2 we introduce some basic 
notation for interest rates and in section 3 we describe the basic requirements 
for LIBOR models. In section 4 we review the construction of LIBOR market 
models, describe its shortcomings and discuss some approximation methods 
developed to overcome them. In section 5 we review forward price models 
and in section 6 we discuss Markov- functional models. Finally, in section 7 
we present afhne LIBOR models and in section 8 we outline the extensions 
of LIBOR models to the multi-currency and default risk settings. 



Let us consider a discrete tenor structure = Tq < Ti < • • • < T/v, with 
constant tenor length 6 = T^+i — T^. The following notation is introduced 
for convenience; K := {1, . . . ,N — 1} and K := {1, . . . , N}. Let us denote 
by B{t,T) the time-t price of a zero coupon bond maturing at time T; by 
L{t,T) the time-t forward LIBOR rate settled at time T and received at 
time T + 6; and by F{t, T, U) the time-t forward price associated to the 
dates T and U. These fundamental quantities are related by the following 
basic equation: 



2. Interest rate markets - notation 



l + 6Lit,T) 



B{t,T) 



F{t,T,T + d). 



(2.1) 



B{t,T + 6) 
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Throughout this work, "B = (Q,J^,F,W) denotes a complete stochastic 
basis, where T = Tt, F = iJ~^t)telo,T]i a-nd T/v < T < oo. We denote by 
AiQP) the class of martingales on "B with respect to the measure P. 

We associate to each date in the tenor structure a forward martingale 
measure, denoted by Pt^.) k £ K. By the definition of forward measures, cf. 
Musiela and Rutkowski (1997, Def. 14.1.1), the bond price with maturity Tk 
is the numeraire for the forward measure Pr^- Thus, we have that forward 
measures are related to each other via 

dPT, 



F{0,Tk,Tk+i) B{0,n) B{t,Tk+i) 
while they are related to the terminal forward measure via 
dPr, 



F{0,Tk,TN) B{0,Tk) B{t,TN)' 



All forward measures are assumed to be equivalent to the measure P. 

3. Axioms for LIBOR models 

In this section, we present and discuss certain requirements that a model 
for LIBOR rates should satisfy. These requirements are motivated both by 
the economic and financial aspects of LIBOR rates, as well as by the practical 
demands for implementing and using a model in practice. The aim here is 
to unify the line of thought in Hunt et al. (2000) and in Keller-Ressel et al. 
(2009). 

A model for LIBOR rates should satisfy the following requirements: 
(Al) LIBOR rates should be non-negative: L{t,T) > for all t G [0, T]. 
(A2) The model should be arbitrage-free: L{-,T) G J^iW^+s)- 
(A3) The model should be analytically tractable, easy to implement and 

quick to calibrate to market data. 
( A4) The model should provide a good calibration to market data of liquid 

derivatives, i.e. caps and swaptions. 
Requirements (Al) and (A2) are logical conditions originating from eco- 
nomics and mathematical finance. Below, we briefly elaborate on (A3) and 
(A4); they stem from practical demands and are more difficult to quan- 
tify precisely. In order to clarify their difference, we point out that e.g. the 
Black-Scholes model obviously satisfies (A3) but not (A4). 

Requirement (A3) means that we can price liquid derivatives, e.g. caps and 
swaptions in "closed form" in the model, so that the model can be calibrated 
to market data in a fast and easy manner. Ideally, of course we would like 
to be able to price as many derivatives as possible in closed form. Here, 
"closed form" is understood in a broad sense meaning e.g. Fourier transform 
methods; really closed form solutions a la Black-Scholes are typically hard to 
achieve. Hunt et al. (2000) say that the model is analytically tractable if it is 
driven by a low- dimensional Markov process. In Keller-Ressel et al. (2009), 
as well as in the present article, we say that a model is analytically tractable 
if the structure of the driving process is preserved under the different forward 
measures. 

Finally, requirement (A4) means that the model is able to describe the 
observed data accurately, without overfitting them. We will not examine this 
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requirement further in this article. On an intuitive level, since the models we 
will describe in the sequel are driven by general Markov processes or general 
semimartingales, we can always find a driving process that provides a good 
calibration to market data. However, an empirical analysis should be per- 
formed in order to identify such a driving process (cf. e.g. Jarrow et al. 2007 
and Skovmand 2008, Ch. III). 



4. LIBOR MARKET MODELS 

LIBOR market models were introduced in the seminal papers of Miltersen 
et al. (1997), Brace et al. (1997) and Jamshidian (1997). In this framework, 
LIBOR rates are modeled as the exponential of a Brownian motion un- 
der their corresponding forward measure, hence they are log-normally dis- 
tributed. This is the so-called log-normal LIBOR market model. Caplets are 
then priced by Black's formula (cf. Black 1976), which is in accordance with 
standard market practice. Later, LIBOR market models were extended to 
accommodate more general driving processes such as Levy processes, sto- 
chastic volatility processes and general semimartingales, in order to describe 
more accurately the market data; cf. Jamshidian (1999), Glasserman and 
Kou (2003), Eberlein and Ozkan (2005), Andersen and Brotherton-Ratcliffe 
(2005), Belomestny and Schoenmakers (2010) and Belomestny et al. (2009), 
to mention just a fraction of the existing literature. 

Consider an initial tenor structure of non-negative LIBOR rates L(0, T^), 
k & K, and let A(-,Tfc) : [0, T^] — )■ M denote the volatility of the forward 
LIBOR rate L[-,Tk), k £ K; the volatilities are assumed deterministic, for 
simplicity. Let H denote a semimartingale on (i7, J^, F, Pt^^) with triplet of 
semimartingale characteristics T(iJ|Pr^) = {B,C,u) and Hq = a.s.; H 
satisfies certain integrability assumptions which are suppressed for brevity 
(e.g. finite exponential moments, absolutely continuous characteristics). The 
process H is driving the dynamics of LIBOR rates and is chosen to have a 
tractable structure under Prjv (s-g- H is a Levy or an affine process). 

In LIBOR market models, forward LIBOR rates are modeled as follows: 

L{t,Tk) = L{0,Tk)exp (^J P{s,Tk)ds + J A(s,Tfc)di/,j , (4.1) 

where (3{-,Tk) is the drift term that makes L{-,Tk) G M{TPTf.^i), for all 
k £ K. Therefore, the model clearly satisfies requirements (Al) and (A2). 
Now, using Theorem 2.18 in Kallsen and Shiryaev (2002), we have that 

/3is,n) = -X{s,Tk)bJ'+' - ^X\s,Tk)cs 

_ J (eA(s,r,)x _ ^ _ x{s,Tk)x)Fl'+\dx), (4.2) 

such that indeed L(-,rfc) G A^(PTfe+J. Here (fef'-'+S c^, fJ'=+') denote the dif- 
ferential characteristics of H under Pt^+i ■ Therefore, in order to completely 
understand the dynamics of the model we have to calculate the characteris- 
tics (6j'=+\c„Fj'=+^). 
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These characteristics fohow readily from Girsanov's theorem for semi- 
martingales (cf. Jacod and Shiryaev 2003, III. 3. 24) once we have the density 
between the measure changes at hand. It is convenient to express this den- 
sity as a stochastic exponential. Keeping (2.3) in mind, and denoting (4.1) 
as follows 

dL{t,Tk) = L{t-,Tk)dH^, (4.3) 

i.e. H'^ is the exponential transform of the exponent in (4.1), we get from 
(2.1) that 



dF{t,n,Tk+i) = 5dL{t,n) 

= SL{t-,n)dH^ 



SL{t-,T, 



l + SL{s-,Tk) 

vO / t 

Therefore, the density between the measure changes takes the form 




5(0, Tfc 



(4.5) 



This calculation reveals the problem of LIBOR market models: the density 
process between the measure changes - and thus the characteristics of H 
under the forward measures - does not depend only on the dynamics of 
H^, or equivalently on the dynamics of f X{s,Ti.)dHs, as is the case in e.g. 
HJM models. It also crucially depends on all subsequent LIBOR rates, as 
the product and the terms in (4.5) clearly indicate. This means, 

in particular, that the structure of the model is not preserved under the 
different forward measures; e.g. if is a Levy or an affine process under 
the terminal measure Ptjvi then H is neither a Levy nor an affine process 
under any other forward measure - not even a time-inhomogeneous version 
of those. Therefore, LIBOR market models do not satisfy requirement (A3). 

The semimartingale H, that drives the dynamics of LIBOR rates, has the 
following canonical decomposition under the terminal martingale measure 

t t 

Ht = Bt + j ^s^Ws + j j x{ii^ - u){ds,dx), (4.6) 

OR 

(cf. Jacod and Shiryaev 2003, II. 2. 38 and Karatzas and Shreve 1991, The- 
orem 3.4.2) where W denotes the P^j^-Brownian motion and //^ denotes 
the random measure of the jumps of H. The P-r^v-compensator of is v 
and C' = /o Cgds. Straightforward calculations using the density in (4.5) 
(cf. e.g. Kluge 2005 or Papapantoleon and Siopacha 2009) yield that the 
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P5"j.^^-Brownian motion VF"^'=+i is related to the Pj^^-Brownian motion via 

L{t-,Ti) 
5L{t-,Ti) 



while the P^^.^^ -compensator of fi^ , , is related to the P^jv -compensator 
of via 

zy^'=+nds,dx) = f n l{s,x,Ti)\ v{dsAx), (4.8) 

\Z=/c+l / 

where 

In addition, the drift term of the LIBOR rate L(-,Tk) relative to the Ptat 
differential characteristics of H, i.e. (J},c,F), is 

V /=fc+i / 

(4.10) 

The consequences of the intractability of LIBOR market models are the 
following. When the driving process is a continuous semimartingale, then 

• caplets can be priced in closed form; 

• swaptions and other multi-LIBOR products cannot be priced in 
closed form; 

• Monte-Carlo simulations are particularly time consuming, since one 
is dealing with coupled high dimensional stochastic processes. 

When the driving process is a general semimartingale, then 

• even caplets cannot be priced in closed form, let alone swaptions or 
other multi-LIBOR derivatives; 

• the Monte-Carlo simulations are equally time consuming. 
Several approximation methods have been developed in the literature in 

order to overcome these problems. We briefly review three of the proposed 
methods below; for other methods and empirical comparison we refer the 
interested reader to the review paper by Joshi and Stacey (2008) . 

4.1. "Frozen drift" approximation. The first and easiest solution to the 
problem is the so-called "frozen drift" approximation, where one replaces the 
random terms in (4.10) or (4.5) by their deterministic initial values, i.e. 

6L{s-,T,) _ 6LiO,Ti) 
l + 5L{s-,Ti) ^ l + 6L{0,Tiy ^ ' ^ 

This approximation was first proposed by Brace et al. (1997), and has been 
used by many authors ever since. Under this approximation, the measure 
change becomes a structure preserving one - observe that the density in (4.5) 
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3^ " 0,15.. 




Figure 1. Difference in implied volatilities between the ac- 
tual LIBOR and the frozen drift prices (left), and between the 
actual LIBOR and the Taylor approximation prices (right). 
Source: Papapantoleon and Siopacha (2009). 

depends now only on the driving process and the volatility structure - and 
the resulting approximate LIBOR market model is analytically tractable; 
e.g. caplets and swaptions can now be priced in closed form even in models 
driven by semimartingales with jumps. 

However, this method yields poor empirical results. Comparing the prices 
of either liquid options, or long-dated options, using the frozen drift approx- 
imation to the prices obtained by the simulation of the actual dynamics 
for the LIBOR rates, we can observe notable differences both in terms of 
prices and in terms of implied volatilities. See Figure 1 for an example. We 
refer to Kurbanmuradov et al. (2002), Siopacha and Teichmann (2010), and 
Papapantoleon and Siopacha (2009) for further numerical illustrations. 

4.2. Log-normal approximations. The following approximation schemes 
for the log-normal LIBOR market model were developed by Kurbanmuradov 
et al. (2002). Consider the log-normal LIBOR market model driven by a one- 
dimensional Brownian motion, for simplicity. The dynamics of LIBOR rates 
(expressed under the terminal measure) take the form 

dL{t,Tk) = L{t,Tk){XtmdWt + f3t{Tk)dt), (4.12) 

where the drift term equals 

A(Tfc) = -A,(r,) JI^^L±IiL_Xtm), (4.13) 

cf. (4.10); w.l.o.g. we can set c = 1. A very crude log- normal approximation 
is to "neglect" the non-Gaussian terms in the SDE, i.e. to set /3f (T^) = 0. Of 
course, this approximation does not yield satisfactory results - in principle, 
results are even worse than the frozen drift approximation. 

One can develop more refined log- normal approximations as follows: let 
/(^) = ' define the process Z, which equals the terms that need to 
be approximated; i.e. 
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Applying Ito's formula, we derive the SDE that Z{-,Tk) satisfies 

dZ{t,Tk) = f'{L{t,n))L{t,n)Xt{Tk)dWt 

+ {f'{L{t,Tk))L{t,Tk)f3t{Tk) + ^f"{L{t,n))L\t,Tk)XUTk)}dt 

=:Akit,Z)dt + Bk{t,Z)dWt, (4.15) 

with the initial condition Z{0,Tk) = /(L(0, T^)). Note that the coefficients 
Ak and can be calculated explicitly, by solving (4.14) for L and substi- 
tuting into (4.15). Moreover, Z in Ak and Bk denotes the dependence on 
the whole vector Z = [Z{-,Ti), . . . , Z{-,T]\f)]. The first and second Picard 
iterations for the solution of this SDE are 

z\t,n) = z{o,n) = (4.16) 

1 + dL{0,Tk) 

and 

t t 
Z\t,n) = Z{0,Tk) + j Ak{s,Z^)ds + j Bk{s,Z'')dWs. (4.17) 



Note that Z^ is constant, while Z^(t,Tk) has a Gaussian distribution since 
the coefficients Ak{-,Z^) and Bk{-,Z^) are deterministic. 

Now, replacing the random terms Z(-,Tfc) in /3{Tk) with the Picard iter- 
ates Z^{-,Tk) and Z'^{-,Tk) leads to two different log-normal approximations 
to the dynamics of LIBOR rates. Obviously, the approximation by Z^ is 
nothing else than the frozen drift approximation. The approximation by Z^ 
is again log-normal, cf. (4.13) and (4.14), and yields very good empirical 
results. This latter approximation has the advantage that the law of the 
approximate LIBOR rate is known at any time t, hence the time-consuming 
Monte Carlo simulations can be avoided. For the empirical and numerical 
analysis of these approximations we refer to Kurbanmuradov et al. (2002) 
and Schoenmakers (2005, Ch. 6). 



4.3. Strong Taylor approximation. Another approximation method has 
been recently developed by Siopacha and Teichmann (2010) and Papapan- 
toleon and Siopacha (2009). The main idea behind this method is to replace 
the random terms in the drift (4.10) by their first-order strong Taylor ap- 
proximations. The Taylor approximation is developed by a perturbation of 
the SDE for the LIBOR rates and a subsequent Taylor expansion. 

Let us denote the log-LIBOR rates by G(-,Tfc) := logL(-,Tfc). Then, they 
satisfy the linear SDE (under the terminal measure) 

dG{t,n)=At,Tk)dt + \{t,n)dHt, (4.18) 

with initial condition G(0, T^) = log L(0, T^); cf. (4.1) and (4.10). We perturb 
this SDE by a real parameter e, i.e. 

dG'{t,Tk) = e0'{t,Tk)dt + \{t,Tk)dHt), (4.19) 

where G^iQ^Tk) = G{0,Tk) for all e. The superscript e in the drift term 
j3^{-,Tk) is a reminder that this term is also perturbed by e, since it contains 
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all subsequent LIBOR rates; see (4.10) again. Now, the first order strong 
Taylor approximation of G*^, denoted by TC^, is 

TG'it,n) = GiO,n) + e^J^^^G'it,n). (4.20) 

We denote the "first variation" process ^|^_qG"^(-, T^.) by y(-,Tfc), and then 
we can deduce, after some calculations, that it has the decomposition 

t t 

Y{t,n) = I ^(s,r,.)ds + J X{s,n)dHs, (4.21) 



where l3^{s,Tk) := /3'^{s,Tk)\e=o- Hence, this is a deterministic drift term, 
obtained by replacing the random terms in (4.10) by their deterministic ini- 
tial values. In particular, we can easily deduce from (4.21) that, for example, 
if i/ is a Levy process then Y{-,Tk) is a time-inhomogeneous Levy process. 
Concluding, we have developed the following approximation scheme: 

logL(t,rfc) « logL{0,Tk) + Y{t,Tk), (4.22) 

where Y{-,Tk) has the decomposition (4.21); compare with (4.1). 

The advantage of this method is threefold: (a) it is universal, and can 
be applied to LIBOR models with stochastic volatility and/or jumps, (b) it 
is faster and easier to simulate than the actual SDE for the LIBOR rates, 
and (c) the empirical performance is very satisfactory; cf. Figure 1 and the 
aforementioned articles for further numerical illustrations. The drawback is 
that it is based on Monte Carlo simulations, hence computational times can 
become long. 



5. Forward price models 

Forward price models were developed by Eberlein and Ozkan (2005), and 
further investigated by Kluge (2005); see also Eberlein and Kluge (2007). We 
consider a setting similar to LIBOR market models, i.e. an initial tenor struc- 
ture of non-negative LIBOR rates, A(-,Tfc) denotes the volatility of the for- 
ward LIBOR rate L(-, T^), and H denotes a semimartingale on ($7, J^, F, IPtat) 
with triplet of characteristics (3,0,1/); again some assumptions on H are 
suppressed. The process H is driving the dynamics of LIBOR rates and has 
a tractable structure under TPtn (^-g- H is a. Levy or an affine process). 

Instead of modeling the forward LIBOR rate directly, one now models the 
forward price in a similar fashion; that is 

l+SL{t,Tk) = (l+5L(0,rfc))exp (^J ^{s,Tk)ds + J A(s,Tfc)d/7,j , (5.1) 

where again the drift term is such that L{-,Tk) G M{^Tu+i)-, for all k £ K; 
i.e. /3{-,Tk) has similar form to (4.2). Therefore the model obviously satisfies 
requirement (A2). 
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Now, from (2.3) and (5.1), we get that the density between the forward 
measures is 



dPT; 



fc + 1 



5(0, Tjv) 




6L{t,Ti)) (5.2) 



^ N-1 







fc+i 



Observe that this density only depends on the driving process H and the 
volatility structures, hence we can deduce that the measure changes between 
forward measures are Esscher transformations; cf. Kallsen and Shiryaev 
(2002) for the Esscher transform. Analogously to eqs. (4.6)-(4.8), we have 
now that the Pr^^j-Brownian motion is related to the P^j^-Brownian mo- 
tion via 

= Wt- f (y V^sds, (5.3) 

^ \i=k+i J 

while the Prj.+i -compensator of fi^ , is related to the P^jy-compensator of 
fi^ via 

/ N-1 \ 
i/^'=+i(ds,dx) = exp X ^ X{s,Ti) ij{ds,dx). (5.4) 
V l=k+l J 

Thus the structure of the driving process is preserved. For example, if is a 
Levy or an affine process under Ptjv> then it becomes a time-inhomogeneous 
Levy or affine process respectively under any forward measure Pt^.^^. This 
implies that requirement (A3) is satisfied, so that caplets and swaptions can 
be priced in closed form. In this class of models, we have the additional 
benefit that we can even price some exotic path-dependent options easily 
using Fourier transform methods; see Kluge and Papapantoleon (2009) for 
an example. 

The main shortcoming of forward price models is that negative LIBOR 
rates can occur, similarly to HJM models, since 

1 + 6L{t, Tfc) > ^ L{t, Tk) > for all t G [0, T^]. 

Therefore, this model can violate requirement (Al). 

6. Markov-functional models 

Markov-functional models were introduced in the seminal paper of Hunt, 
Kennedy, and Pelsser (2000). In contrast to the other approaches described 
in this review, the aim of Markov-functional models is not to model some 
fundamental quantity, e.g. LIBOR or swap rates, directly. Instead, Markov- 
functional models are constructed by inferring the model dynamics, as well 
as their functional forms, through matching the model prices to the market 
prices of certain liquid derivatives. That is, they are implied interest rate 
models, and should be thought of in a fashion similar to local volatility 
models and implied trees in equity markets. 
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The main idea behind Markov-functional models is that bond prices and 
the numeraire are, at any point in time, a function of a low- dimensional 
Markov process under some martingale measure. The functional form for 
the bond prices is selected such that the model accurately calibrates to the 
relevant market prices, while the freedom to choose the Markov process 
makes the model realistic and tractable. Moreover, the functional form for 
the numeraire can be used to reproduce the marginal laws of swap rates or 
other relevant instruments for the calibration. 

More specifically, let (M, M) denote a numeraire pair, and consider a 
(time-inhomogeneous) Markov process X = {Xt)o<t<T under the measure 
M. In the framework of Markov-functional models, one assumes that bond 
prices have the functional form 

B{t,S) = B{t,S;Xt), 0<t<ds<S, (6.1) 

where ds denotes some "boundary curve". In applications, the boundary 
curve typically has the form 

where T^, is a common time of maturity. One further assumes that the nu- 
meraire M is also a function of the driving Markov process X, i.e. 

Mt = M{t; Xt), 0<t<T. (6.3) 

Therefore, in order to specify a Markov-functional model, the following quan- 
tities are required: 

(PI) the law of X under the measure M; 

(P2) the functional form B{ds, S; •) for S G [0, T]; 

(P3) the functional form M{t; •) for < t < T. 

In applications, the Markov process is specified first and is typically a 
diffusion process with time-dependent volatility. Then, the functional forms 
for the bond prices and the numeraire are implied by calibrating the model 
to market prices of liquid options. The choice of the calibrating instruments 
depends on the exotic derivative that should be priced or hedged with the 
model. If the exotic depends on LIBOR rates, e.g. the flexible cap, then the 
model is calibrated to digital caplets, which leads to the Markov- functional 
LIBOR model. If the exotic depends on swap rates, e.g. the Bermudan swap- 
tion, then the model is calibrated to digital swaptions, which leads to the 
Markov-functional swap rate model. Let us point out that the functional 
forms are typically not known in closed form, and should be computed nu- 
merically. These models typically satisfy requirements (Al), (A2) and (A3). 
For further details and concrete applications we refer the reader to the books 
by Hunt and Kennedy (2004) and Fries (2007), and the references therein. 

Remark 6.1. Let us point out that forward price models and affine LIBOR 
models, that will be introduced in section 7, belong to the class of Markov- 
functional models, while LIBOR market models do not. In LIBOR market 
models the LIBOR rates are functions of a high- dimensional Markov process. 
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6.1. Markov- functional LIBOR model. In order to gain a better un- 
derstanding of the construction of Markov-functional models, we will briefly 
describe a Markov-functional model calibrated to LIBOR rates. This model 
is called the Markov- functional LIBOR model. 

The set of relevant market rates are LIBOR rates L{-,Tk), k £ K. We will 
consider the numeraire pair (M, M) = T/v), P^-^). 

In order to be consistent with Black's formula for caplets, we assume that 
L{-,Tn-i) is a log-normal martingale under Ptad i-G- 



dLit,TN- 



ait,TN^i)L{t,TN^i)dWt. 



(6.4) 



where W denotes a standard Brownian motion under Ptjv (t(-,T/v-i) 
is a deterministic, time-dependent volatility function. Hence, we have that 



L{t,TN-i) = L(0,rAr„i)exp 



-St + Xti, 



(6.5) 



where S = Jq cr^(s, T/v-i)ds, and X is a deterministic time-change of the 
Brownian motion, that satisfies 



dXt = a{t,TM-i)dWt. 



(6.6) 



We will use X as the driving process of the model, which specifies (PI). 
Regarding (P2), the boundary curve is exactly of the form (6.2) with 
= Tjv-i, hence we need to specify B{Ti,Ti; XtJ for i £ K, which is 
trivially the unit map. We also need to specify B(T/v-i, T^r; J; using 
(2.1) and (6.5) we get that 

Then, we can recover bond prices in the interior of the region bounded by 
ds using the martingale property: 



B{t,S;Xt) = B{t,TN;Xt)JETj, 



B{ds,S;X9^] 



6(55, Tat; 



{61 



Now, it remains to specify the functional form B(Tj, T/v; X^-), i £ K, for 
the numeraire, cf. (P3). In the framework of the Markov- functional LIBOR 
model, this is done by deriving the numeraire from LIBOR rates and in- 
ferring the functional forms of the LIBOR rates via calibration to market 
prices. Equation (2.1) combined with (6.2) and the fact that i?(Tj,Tj+i) is a 
function of Xt^, cf. (6.8), yield that L{Ti,Ti) is also a function of Xjv . The 
functional form is 

l + 6L{T,,Ti;XT,) ^ 



^iTi,Ti+i;XTi] 



1 



B{Ti,TN; XtJ'S^t^ 



_B(Ti+i,Tjv;XTj^-^) 

Rearranging, we get the following functional form for the numeraire 

1 



B(Tj, Tat; X^-J 



il + 6L{Ti,Ti;XTj)JET^ 



(6.9) 



(6.10) 
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This formula provides a backward induction scheme to calculate B(rj, T/v; •) 
from B(Tj_|_i, T/v; •) for any value of the Markov process; the induction starts 
from B{Tn,Tn) = 1. 

The calibrating instruments are digital caplets with payoff ^{L{Ti,Ti)>K}i 
i E K, and their market values are provided by Black's formula; we de- 
note them by Vo(Tj,/C). Assuming that the map ^ i— )■ L(Tj,Tj;^) is strictly 
increasing, there exists a unique strike lC{Ti,x*) such that the set equality 

{Xt, > X*} = {L{Ti,Ti;XT,) > }C{Ti,x*)} (6.11) 

holds almost surely. Define the model prices 



,m,x*) = B{0,TN)ETr. 



-1{Xt,>x*} 



(6.12) 



B(Tj, T/v; Xxi 

which have to be calculated numerically. Therefore, we can equate market 
and model prices 

Yom,IC{Ti,x*))=Vo{Ti,x*), (6.13) 

where the strike IC{Ti,x*) is determined by Black's formula using some nu- 
merical algorithm. 

Hence, we have specified, numerically at least, the functional form for the 
LIBOR rates, which yields also the functional form for the numeraire via 
(6.10). This completes the specification of the Markov-functional LIBOR 
model. This model satisfies requirement (A3), in the sense of Hunt et al. 
(2000), since all bond prices are functions of a one-dimensional diffusion. 

7. Affine LIBOR models 

Affine LIBOR models were recently developed by Keller-Ressel, Papa- 
pantoleon, and Teichmann (2009) with the aim of combining the advantages 
of LIBOR market models and forward price models, while circumventing 
their drawbacks. We provide here a more general outline of this framework, 
which is based on two key ingredients: martingales greater than 1, which are 
increasing in some parameter. 

The construction of martingales greater than 1 is done as follows: let 
be an J^^^-measurable, integrable random variable, taking values in [l,oo), 
and set 

Mj" = E[y^|Ji], 0<t<T. (7.1) 

Then, using the tower property of conditional expectations, it easily follows 
that M'^ = (M")o<f<T is a martingale. Moreover, it obviously holds that 
> 1 for ah t G [0,T]. 
In addition, assume that the map u i— )• Yji is increasing; then, we imme- 
diately get that the map 

u ^ (7.2) 

is also increasing, for all t G [0, T]. 

Now, using the family of martingales we can model quotients of bond 
prices as follows. Consider a decreasing sequence {uk)f^^-j^ and set 



B{t,TN) * 



M^", tG[0,Tk],kGK, (7.3) 
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requiring that the initial values of the martingales fit today's observed mar- 
ket prices, i.e. q^q''^'^-^ = Mq'' . Since is increasing in u, we have that 

M^" > M^' for k<l^ Uk > ui. (7.4) 

Hence, we can deduce that bond prices are decreasing as functions of time 
of maturity, i.e. B(t,Tk) > B{t,Ti) for k <l. 

Turning our attention to LIBOR rates, we get that 



B{t,Tk) _ M^" 



1 + 6Lit, n) = Tjyf^ = > 1' (7-5) 

for all t £ [0, Tfc] and all k £ K; this is a trivial consequence of (7.4). More- 
over, the martingale property of the LIBOR rate under its corresponding 
forward measure follows easily from the structure of the measure changes 
(2.3), and the structure of the martingales. Indeed, we have that 

M^k+i dPT^ M"fc+i Mo'=+' ^ ^ ' 

Therefore, we have just described a broad framework for modeling LIBOR 
rates, in which requirements (Al) and (A2) are satisfied. The next step is to 
show that requirement (A3) is also satisfied. We will not pursue this in full 
generality, instead we will consider a specific form for the variable 1^, and 
thus for the martingales M". In addition, the model is driven by an affine 
process, and is henceforth called the affine LIBOR model. 

7.1. AfRne processes. Let X = {Xt)o<t<T be a stochastically continuous, 
time-homogeneous Markov process with state space D = M^q, starting from 
X £ D. The process X is called affine if the moment generating function 
satisfies 

E4e<"'^*>] = exp {Mu) + {Mu),x)), (7.7) 

for some functions </> : [0, T] x M and : [0, T] x Xt K'', and all 

{t,u,x) £ [0, T] xIt y. D, where 

It := G M'^ : [e<"'^^>] < 00, for all x £ (7.8) 

We will assume in the sequel that £ I^. The functions (p a-nd ■0 satisfy the 
semi-flow property 

4>t+s{u) = Mu) + (psiM'^)) , 

(7.9) 

with initial condition 

(j)o{u) = and ipo{u) = u, (7-10) 

for all {t+s,u) £ [0, T] xXy. Equivalently, (j) and ^p satisfy generalized Riccati 
differential equations. For comprehensive expositions of affine processes we 
refer the reader to Duffie et al. (2003) and Keller-Ressel (2008). 
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7.2. Affine LIBOR model. In the affine LIBOR model, the random vari- 
able Yji has the following form: 

YU^^{u,Xt)^ (7.11) 

where u G IR>q and Xt is a random variable from an M>Q-valued afhne 
process X. Hence, > 1, while the map u i— )• Yji is obviously increasing; 
note that inequalities involving vectors are understood componentwise. 

Using the Markov property of affine processes, we can deduce that the 
martingales have the form 

Mj" = ]E[e<"'^^>|J-j] 

= exp {^T-tiu) + {^PT-tiu),Xt)). (7.12) 

Therefore, LIBOR rates have the following evolution: 

l + 5L{t,n) = j^=exp{Ak,t + {Bk,t,Xt)), (7.13) 



where 



^k,t ■= 4>TN-t{uk) - 4>Tj^-t{uk+l) 
Bk,t ■■= IpTN-tiUk) - tpTN-t{uk+l). 



(7.14) 



Let us point that, under reasonable assumptions on the driving affine pro- 
cess, we can prove that the affine LIBOR model can fit any term structure 
of initial LIBOR rates; cf. Proposition 6.1 in Keller-Ressel et al. (2009). 

Now, regarding requirement (A3), let us turn our attention to the struc- 
ture of the driving process under the different forward measures. Using the 
connections between forward measures (2.3), the Markov property of affine 
processes, and the flow equations (7.9), we can show that 

Et, [e^-'^'-^lJ-,] = Et, [^e^"'^") | J",] (7.15) 



exp (^(pr-silpTN-riuk) + v) - (j)r-silpTN-r{uk)) 

+ (ipr-silpTN-riUk) + v) - (V'Tjv-r (^ifc)) , A^) ^ : 



cf. Keller-Ressel et al. (2009, eq. (6.15)). This means that X becomes a time- 
inhomogeneous affine process under any forward measure. Note that the 
measure changes are again Esscher transformations, similarly to forward 
price models. Consequently the affine LIBOR model satisfies requirements 
(Al), (A2) and (A3). 

The pricing of caplets and swaptions in the affine LIBOR model using 
Fourier transform methods is described in Keller-Ressel et al. (2009). In 
addition, closed-form valuation formulas ~ in terms of the ^^-distribution 
function - are derived when the driving affine process is the Cox~Ingersoll- 
Ross (CIR) process. 

8. Extensions 



The different approaches for modeling LIBOR rates have been extended 
in two different directions: (i) to model simultaneously LIBOR rates for 
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different currencies and the corresponding foreign exchange rates, and (ii) 
to jointly model default-free and defaultable LIBOR rates. 

8.1. Multiple currencies. The log-normal LIBOR market model has been 
extended to a multi-currency setting by Schlogl (2002) and by Mikkelsen 
(2002). The Levy LIBOR model and the Levy forward price model have 
been extended to model multiple currencies and foreign exchange rates by 
Eberlein and Koval (2006) . A multi-factor approach to multiple currency LI- 
BOR models has been presented in Benner et al. (2009). Markov-functional 
models have been extended to the multi-currency setting by Fries and Rott 
(2004) and Fries and Eckstiidt (2010). 

8.2. Default risk. The log-normal LIBOR market model has been first ex- 
tended to model default risk by Lotz and Schlogl (2000), who used a deter- 
ministic hazard rate to model the time of default. Eberlein et al. (2006), 
borrowing also ideas from Schonbucher (2000), constructed a model for 
default-free and defaultable rates where they use time-inhomogeneous Levy 
processes as the driving motion and the "Cox construction" to model the 
time of default (cf. e.g. Bielecki and Rutkowski (2002) for the Cox construc- 
tion). This has been extended to a model where defaultable bonds can have 
rating migrations by Grbac (2010). 
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